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Based on the branching fractions of J/t[>(ip(2S)) — > VP from different collaborations, the pseudoscalar mixing 
is extensively discussed with a well established phenomenological model. The mixing angle is determined to be 
— 14° by fitting to the new world average if only quark content is considered. After taking into account the 
gluonic content in r\ and r\ simultaneously, the investigation shows that rj favors only consisting of light quarks, 
while the gluonic content of 77' is Z^, = 0.30 ± 0.24. 

PACS numbers: 13.25.Gv,14.40.Ag 

1 Introduction 

As the ground pseudoscalar nonet, n, K, 77 and 77', in the constituent quark model, their masses and widths 
are determined with high precision and the main decay modes are also observed[I] in addition to the forbidden 
and rare decays. However there is one issue, pseudoscalar mixing, remains not completely settled, which has been 
discussed for many times with different transitions. The linear Gell-Mann-Okubo(GMO) mass relation[2] gives a 
mixing angle, dp = —11°, which is hardly consistent with the value, Bp = —24.6°, obtained from the quadratic 
GMO mass formula by replacing the meson masses by their squares. The full set of J/ip decays into a vector and 
a pseudoscalar was measured by Marklll, and the phenomenological analysis of mixing angle is determined to be 
9p = (—19.2 ± 1.1: 3 . which was confirmed by DM2 1 . Both of them got the conclusion that 77 and 77' consist 
of light quarks, with no contribution from gluonium or radial excitation states. After that an important work was 
performed by Bramon and Scadron 5, 6 , taking into account w-<p mixing in the analysis for J/ip —¥ VP, a weighted 
dp is calculated to be (—15.5 ± 1.3)° based on many different transitions. For a nice review based on the discussions 
before 2000, see[7] in which the reasonable range of 77-77' mixing angle is believed to be —20° ~ —10°. 

Recently the new experimental data on J/ip -4 VP and ip(2S) -> VP were reported by BES P 19] FTO l ITT ] [12 ^ [T3 l 
[14], BABAR[15I [TBI [T7l H~8] and CLEOjTS]. It is worth pointing out that some of the new measurements are not 
well consistent with the previous works. Take J/ip — > pir for example, the branching fractions measured by BES is 
(2.10±0.12)% 8 , subsequently confirmed by BABAR|15], which is larger than the world average (1.28±0.10)%[20], 
about 64%. This significant change stimulates new interest in this issue [21] [22] [23] [24]. The analysis in Ref. [23] 
indicates it is difficult to get reasonable results with the updated branching fractions of J/ip — > pn, however, the 
results in Ref. [24] performed with the same data and phenomenological model seems reasonable. This discrepancy 
motivated us to reanalyze the full set of J/tp — > VP data. Actually it is difficult for us to compare the results 
obtained with different sets of parameters in one time, in this paper we would like to discuss this issue for different 
cases, e.g. fix SU(3) breaking term x to 0.64, 0.82 or 1. 
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2 Notation 



The physical eigenstates 77, rf are the mixture of octet, singlet and gluonium. And they are defined as, 

\n >= X V \N > +Y n \S > +Z n \G >, 
W >= X' V \N > +Y^\S > +Z' n \G > . 

where, N — -^={uu + dd), S — ss and G for gluonium; Xi, Yi and Zi denote the magnitude of non-strange, strange 
contents and gluonium in 77 and rf . 

The above form can be written in terms of the three Euler angles, with 

X rj = cos 4> p cos 4>ci ■ 
Y v = -sin0 p cos0Gi, 

Z n = -sin^Gi, ^ 
Xrji = cos cj>p cos4>G2 — siiiipp sin <pG2 sin^ci) 
Y n > = sin 4>p cos 4>G2 + cos 4>p sm.4>G2 sinc^ci, 
Z n > = -sin0G2Cos</> G i. 

If we only consider the simplest case and neglect possible mixing of the rj and 77' with other pseudoscalar states, 
77-77' mixing is characterized by a single mixing angle dp. 

\r\ >= cos9p\i] 8 > -sm9 P \i] >, 

V J 

1 77 >= sin^p^s > + cos 9 P \r/ > . 

where rj and 77' are the orthogonal mixture of the respective singlet and octet iso-spin zero states. 770 and 77s are 
SU(3) quark basis states which are denoted as 770 = "^Iww + dd + ss > and r/g = -^=\uu + dd — 2ss > respectively. 

In terms of quark basis, the 77 and 77' include non-strange and strange contents. In the flavor SU(3) quark model, 
they are defined through quark-antiquark(<7<7~) basis states as, 



X v = Y v > = y ^ cos 9 P - d ^ sin 9 p = cos tf> P , 
X v > = -Y n = sin 9 P + Jl cos9 P = sm<j) P . 



(4) 



where, 9 P = <p P - 54.7°. 



3 Phenomenological model 

J/ip and ip(2S) have the similar decay mechanism and are suppressed by Okubo-Zweig-Iizuka (OZI) rule. Both 
of them decay into a vector and pseudoscalar meson via three gluon annihilation and electromagnetic decays. 
Therefore, in this paper, the phenomenological model for J/tp — > VP in Ref.[25] is simply applied in ip(2S) decays 
to discuss the 77-77' mixing and other physics. 

A first-order parameterization of the amplitudes appears in Ref.|25] and is described there in detail. The 
amplitude, which has contributions from both the three gluon annihilation and electromagnetic processes can be 
expressed in terms of an SU(3) symmetric single-OZI(SOZI) amplitude <?, an electromagnetic amplitude e (the 
coupling strength e has a relative phase 9 e to the strength g because these are produced from different origins) and 
the nonet-symmetry-breaking double-OZI(DOZI) amplitude r, relative to g. SU(3) violation has been accounted for 
by a pure octet SU(3) breaking term. The SU(3) breaking term in strong interaction and electromagnetic process 
are expressed by (1-s) and x, respectively. A factor (1-s) for every strange quark contributing to g and a factor for 
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x for a strange quark contributing to e. The factor s v (s p ) is for the strange vector (pseudoscalar) contributing to r. 

In spite of these simplified assumptions this phenological model contains a rather large number of parameters 
(<7, e, r, s, s p , s v , x, Op and e ). This x can be well determined via V — > P7 and P — > T^7 data, we reanalyzed 
it using the phenomenological model in Ref.[35] and the branching fractions of V — > P7 and P — » 1^7 in Ref.Q], 
x is determined to be 0.82 ± 0.05 and 0y = (3.2 ± 0.9)° and Op = (—12.9 ± 0.5)° which are in good agreement 
in those in Ref.|21). To further simplify it again, s p is ignored in this paper and s v is discussed below with two 
assumptions (s v — and s v = s). 



Table 1: Branching fractions of J/ip -t VP(xl0" 4 ) 



Decay Modes 


Marklll 




DM2 


BES 


BABAR PDG2010 


pir 


142 ± 1 ± 9 




132 ±20 


210 ± 12 ±20.1 


218 ±19 169 ±15 


pn 


1.93 ± 0.13 ±0. 


.29 


1.94 ±0.17 ±0.29 




1.93 ±0.23 


pn' 


1.14 ± 0.14 ±0. 


.16 


0.83 ±0.30 ±0.12 




1.05 ±0.18 


<f>7T 


< 0.068 






< 0.064 


< 0.064 




6.61 ± 0.45 ±0. 


.78 


6.4±0.4± 1.1 


8.98 ±0.24 ±0.89 


14±6±1 7.5 ±0.8 


H 


3.08 ± 0.34 ±0. 


.36 


4.1 ±0.3 ±0.8 


5.46 ±0.31 ±0.56 


4.0 ±0.7 


con 


4.82 ± 0.19 ±0. 


.64 


3.6 ±0.28 ±0.54 


5.38 ±0.12 ±0.65 


4.5 ±0.5 


uir] 


17.1 ±0.8 ±2. 


.0 


14.3 ± 1.0 ±2.1 


23.52 ±2.73 


14.4 ±4.0 ±1.4 17.4 ±2.0 


LOT]' 


1.66 ± 0.17 ±0. 


.19 


l-8±o.8±0.3 


2.26 ±0.43 


1.82 ±0.21 


K*-K+ + c.c. 


52.6 ± 1.3 ±5. 


.3 


45.7 ±1.7 ±7.0 




52±4±1 51.2 ±3.0 


K*°K° + c.c. 


43.3 ±1.2 ±4. 


.5 


39.6 ± 1.5 ±6.0 




48±5±1 43.9 ±3.1 



4 Results 

The experimental data sets shown in Table Q] are analyzed with the least squares method to determine the 
coupling strengths and mixing angle. To clarify the results obtained from different data set, we divided it into 
several subsections to investigate the pseudoscalar mixing. 

4.1 Analysis of J/ip VP from Marklll and DM2 

We start to perform the fit to experimental data with the simplest case, the oj-(p mixing and gluon content are 
ignored. Actually the treatment on SU(3)-breaking parameter x and the second order corrections s v in Ref.[3] and 
Ref.[4] are different. The x is set to 1 and the correction terms s v — 0, is ignored in Marklll analysis, while x 
is fixed to 0.64 and the correction terms s v — s is included in DM2 analysis. To clearly compare the difference 
between them, all the possible combinations are considered to perform the fit. 

A fit to the data without considering SU(3) breaking as well as in Marklll analysis yields 9p = (—13.95 ± 2.39)° 
with x 2 /d-°-f — 9.0/4, which is obviously inconsistent with the value (— 19. 2± 1.4)° [5]. After tuning the parameter, 
we also get a reasonable results which are the same as those in Ref.[5], 

g = 1.10 ± 0.03, s = 0.12 ± 0.03, e = 0.122 ± 0.005, e = 1.25 ± 0.12, p = (-19.34 ± 1.40)°, r = -0.15 ± 0.09. 
But the goodness of fit, x 2 jd.o.f = 10.1/4, seems slightly worse. Compared with the results listed in the first 
column of Tabled s and r also change significantly. The results of the fits performed with x = 0.64 and x = 0.82 
are also given in Tabled Apart from the mixing angle, the values of other parameters are also consistent with the 
previous fit. 

If s v is replaced with s and x is fixed to 0.64, the fit gives Op = (—18.59 ± 1.40)° with x 2 /d.o.f = 9.0/4 which 
is in good agreement with DM2's result Op — (—19.1 ± 1.4)°. Meanwhile we also checked the fits with x = 1 and 
x = 0.82 and the results are listed in Table[2] Compared with the results without considering the contribution of s v , 
the results change significantly, in particular for s, r and p . This is reasonable because the two phenomenological 
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models are slightly different. The fit to DM2 data is also performed to check the discrepancy discussed above. 
In the DM2's analysis, the common error of the branching fractions is removed, so the fitting error here is larger 
than those in Ref.|4j. Here it is clear that the reasonable results can also be obtained Op = (—14.84 ± 4.35)°, with 
X 2 /d.o.f = 1.9/4 in the case of s v = and x = 1. 

Based on the above results, we can get the conclusion that s v plays an important role in the fit to extract the 
mixing angle. The mixing angle in DM2 analysis is consistent with that in Marklll because the latter is not from 
the best fit. 

4.2 Analysis of J/ijj VP from BES, BABAR and PDG2010 

Until now the pseudoscalar mixing is investigated with the well established models and the data measured about 
20 years ago. The new measurements reported by BES, BABAR and the new world average of 2010 are listed in 
TableHJ Each branching fraction is regarded as one constraint in the fit to BES and BABAR data. The amplitude 
of J/tp — > pr\ and J/ip — > pr)' is removed from the fit because no new measurements are available. The results of 
the fits with s v — yields the mixing angle P ~ —17°, which is still consistent with the above results within one 
standard deviation. This value is also in agreement with the previous work in Ref. [2TJ [22J [24] . The fit with s v = s 
is performed, but the quality of fit is very poor. 

The further check is performed using the world average of 2010 [I], and the results are shown in Table [3] As we 
expected, the results are fine for the fit with s v = and the mixing angle p favors ~ 14° . The goodness of the fit 
with s v = s is still worse because of the weight of new measurements in the world average. 

4.3 Analysis of ijj{2S) -> VP 

We now turn to examine the full set of "0(25') — > VP to get the pseudoscalar mixing using the same phenomeno- 
logical model. At present the measurements of i^(2S) — s> VP mainly come from BES and CLEO's reports which 
are shown in Table 2] We have omitted the known upper limit for the i/i(2S) — > (fm and 4>(2S) — > ujt] branching 
fractions in our analysis because they are the upper limits at 90% confidence level rather than branching fractions. 
As previously stated, we just consider the mixing angle between r\ and rf and assume the mixing of lo and <fi is ideal. 
The results listed in Table [5] indicate that both of the above two slightly different models are reasonable and the P 
favors ~ —12° with large uncertainty. Without considering the branching fraction of ip(2S) — > pir, the fit was also 
performed in Ref. |24j . the mixing angle is calculated to be — 10^g which is in agreement with our result. But the 
branching fraction of ip(2S) — > K* + K~, (8.5 ± 4.0) x 10~ 5 , applied in the analysis is not correct. Therefore the 
values of parameters listed in Table [S] are inconsistent with those in Ref. [23] . 







Table 2: 


Results of fit to Marklll data 






Parameter 


s v = 0,x=l 


s v = 0,x=0.64 


s v = 0,x=0.82 


s v = s,x=l 


s v = s,x=0.64 


s v = s,x=0.82 


g 


1.30 ±0.04 


1.31 ±0.04 


1.30 ±0.04 


1.12 ±0.04 


1.11 ±0.04 


1.11 ±0.04 


s 


0.27 ±0.02 


0.28 ±0.02 


0.27 ±0.02 


0.13 ±0.03 


0.13 ±0.02 


0.13 ±0.03 


e 


0.124 ±0.005 


0.123 ±0.05 


0.124 ±0.05 


0.123 ±0.005 


0.123 ±0.005 


0.123 ±0.005 


e 


1.21 ±0.12 


1.29 ±0.12 


1.27±0.12 


1.27 ±0.12 


1.30 ±0.12 


1.29 ±0.12 


Y 


-0.37 ±0.01 


-0.37 ±0.01 


-0.37 ±0.01 


-0.16 ±0.01 


-0.15 ±0.01 


-0.15 ±0.01 


0p 


-13.95 ±2.39 


-13.17 ±2.40 


-13.49 ±2.38 


-18.29 ± 1.43 


-18.59 ± 1.40 


-18.47 ± 1.41 


X 2 /d.o.f 


9.0/4 


7.9/4 


8.3/4 


8.1/4 


9.0/4 


8.6/4 



4.4 uj-(f) mixing 

In the above analysis, the uj-cj) mixing is ignored to simplify the model. This fit in the case of s v = and x — 0.82 
is an attempt to account for the uj-<f> mixing. If u!-<fi mixing angle is left as a free parameter, the fit to the world 
average of 2010 leads to a minimum x 2 = 3.3 for three degrees of freedom, 
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Table 3: Results of fit PDG2010 data 



Parameter 


s v = 0,x=f 


Sv = 0,x=0.64 


Sv = 0,x=0.82 


Sy - — S 1 


s v = s,x=0.64 


s v = s,x=0.82 


g 


1.35 ±0.04 


1.36 ±0.04 


1.36 ±0.04 


1.15 ±0.04 


1.14 ±0.04 


1.14 ±0.04 


s 


0.30 ±0.02 


0.30 ±0.03 


0.30 ±0.02 


0.15 ±0.03 


0.14 ±0.03 


0.14 ±0.03 


e 


0.120 ±0.005 


0.121 ±0.04 


0.121 ±0.04 


0.119 ±0.005 


0.119 ±0.005 


0.119 ±0.005 




1.31 ±0.12 


1.36 ±0.12 


1.34 ±0.12 


1.35 ±0.12 


1.38 ±0.12 


1.36 ±0.12 


r 


-0.37 ±0.01 


-0.37 ±0.01 


-0.37 ±0.01 


-0.16 ±0.01 


-0.15 ±0.01 


-0.15 ±0.01 


Op 


-14.27 ±2.44 


-13.90 ±2.35 


-14.04 ±2.37 


-17.66 ± 1.81 


-17.96 ± 1.77 


-17.84 ± 1.78 


X 2 /d.o.f 


3.1/4 


3.5/4 


3.3/4 


16.5/4 


18.1/4 


17.4/4 



Table 4: Branching fractions of ip(2S) -> VP(xlO~ 5 ) 



Decay modes 



BES 



CLEO 



PDG2010 



fm 
PV 

pi 



)7T 



UJTT 
0J1] 

uirj' 

K*-K+ + c.c. 
K*°K° + c.c. 



5.1 ± 0.7± 1.1 2.4 ±0.8 ±0.2 



1.78t°;g±0.17 
l.87±\H ± 0.33 

< 0.4 

3.3 ± 1.1 ±0.5 
3.1 ± 1.4±0.7 
1.87±^ ±0.28 

< 3.1 
3.2^ ±0.7 
2.9t^±0.4 
13.3±1;| ± 1.7 



3.0± O g ±0.2 



< 0.7 

2.0±i;f =b 0.4 

2.5±i o ± 0.2 

< 1.1 



1.3±J;? ± 0.3 
9.2±i 7 2 ±0.9 



3.2 ± 1.2 
2.2 ±0.6 



+ 1.7 
-1.2 



1.9 

< 0.4 

9 S+ 1 ' 

3.1 ±1.6 
2.1 ±0.6 

< 1.1 

3.2±l:« 

i 7+O.8 
i -'-0.7 

10.9 ±2.0 



Table 5: Results of fit to PDG2010 data of ijj(2S) -> VP 



Parameter 


s v = 0,x=l 


Sv = 0,x=0.64 


s v = 0,x=0.82 


Sy - — tS^X 1 


s v = s,x=0.64 


s v = s,x=0.82 


g 


0.64 ±0.11 


0.65 ±0.10 


0.65 ±0.10 


0.64 ±0.11 


0.65 ±0.04 


0.65 ±0.04 


s 


0.003 ±0.18 


-0.01 ±0.19 


-0.05 ±0.18 


0.02 ±0.18 


-0.10 ±0.19 


-0.10 ±0.20 


e 


0.23 ±0.02 


0.23 ±0.02 


0.23 ±0.02 


0.23 ±0.02 


0.23 ±0.02 


0.23 ±0.02 


0e 


2.73 ±0.62 


2.81 ±0.60 


2.79 ±0.63 


2.75 ±0.64 


2.83 ±0.64 


2.83 ±0.62 


r 


0.18 ±0.28 


0.17 ±0.27 


0.16 ±0.28 


0.14 ±0.28 


0.14 ±0.29 


0.14 ±0.31 


Op 


-12.07 ± 10.42 


-11.94 ± 10.48 


-11.99 ± 10.46 


-11.80 ±10.63 


-12.19 ± 11.59 


-12.19 ±12.18 


xVd.o.f 


4.4/3 


4.5/3 


4.4/3 


4.4/3 


4.5/3 


4.5/3 
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g = 1.36 ± 0.04, s = 0.30 ± 0.03, e = 0.121 ± 0.005, e = 1.33 ± 0.12, 9 P = (-14.06 ± 2.37)°, r = -0.37 ± 0.02, 
6 V = (0.09 ±4.13)°. 

If we assumed s v = s, then the fit with \ 2 /d.o.f of 17.4/3 gives, 

g = 1.14 ± 0.05, s = 0.14 ± 0.04, e = 0.119 ± 0.005, 9 e = 1.36 ± 0.12, 9 p = (-17.78 ± 2.70)°, r = -0.15 ± 0.01, 
9 V = (0.11 ±3.78)°. 

9v is very close to zero and the uncertainty is very large compared with other parameters. This means that there 
is not a significant constraint on it. Among J/tp — »■ VP decays, the amplitude of J/tjj — > <pTr° is directly related to 
the u-(f> mixing, but it is still not observed yet. No observation of J/tp — > (f>ir shows that the contribution of ur-cj) 
is small. On the other hand, the values of other parameters are almost the same as those listed in Table [3] without 
considering uj-(f> mixing. Therefore it is reasonable that u-(p mixing is assumed to be ideal and could be ignored 
in the above analysis. Further check is done by fixing the u> — <f> mixing angle to 3.2° obtained from V — > jP and 
P —> -fV process. The fit with x 2 /d.o.f = 3.3/4 gives, 

g = 1.36±0.04, s = 0.30±0.02, e = 0.121±0.004, 9 e = 1.33±0.13, 9 P = (-14.05 ± 2.36)°, r = -0.37±0.01, these 
values are also in good agreement with those in the hypothesis of the ideal oj-<f> mixing. 

4.5 Gluon content in 77 and rf 

At present 77 is believed to be well-understood as an SU(3) flavor octet with a small quarkonium singlet admixture, 
and not much room for a significant gluonium admixture [3TJ [53] ■ Therefore the analyses [23] are usually performed 
to determine the gluonic content in 77' with the assumption of no gluonic content in 77. After taking into account 
the gluonic content in 77 and 77' simultaneously, we present the fit with the above two slightly different models. In 
the first CclSG, Sy IS assumed to be zero and the fit to the world average in 2010 yields, 

g = 1.32 ± 0.06, s = 0.27 ± 0.04, e = 0.126 ± 0.007, 9 e = 1.34 ± 0.12, 9 P = (-10.21 ± 4.48)°, r = -0.45 ± 0.08, 
4> g i = 0.04 ± 0.05, <j) g2 = 0.53 ± 0.24, r' = -0.77 ± 0.46, 
with x 2 /d.o.f = 1.56/1. 

The second fit is performed under the hypothesis of s v — s, and the results with \ 2 /d.o.f = 3.5/1 are listed as 
follows, 

g = 1.28 ± 0.06, s = 0.24 ± 0.03, e = 0.128 ± 0.007, 9 e = 1.35 ± 0.11, 9 p = (-9.17 ± 4.67)°, r = -0.67 ± 0.08, 
4> a i = 0.11 ± 0.04, g2 = 0.50 ± 0.22, r' = -0.85 ± 0.56. 

The goodness of the second fit is still worse than the first fit. Based on the results of the first fit, the magnitudes 
of gluon components in 77 and 77' are calculated to be = 0.002 ± 0.002 and Z"i, = 0.30 ± 0.24, respectively. The 
small gluonic contribution in 77 shows there is not much room for gluonium admixture, which is consistent with 
the results presented in Ref . [2"2"] . It seems that 30% of 77' component could be attributed to gluonium, but further 
investigation with more precisely data needed to be done due to the large uncertainty. 

5 Summary and outlook 

A wide set of data on J/tp — > VP and tp(2S) — > VP decays are analyzed in terms of a rather general phe- 
nomenological model in an attempt to determine the magnitudes of components in 77 and 77'. The data include the 
branching fractions of J/ip — > VP which were measured nearly 20 years ago and the recent measurements by BES 
and BABAR. The measurements of Marklll and DM2 are reanalyzed, we found that the results obtained from the 
two different phenomenological models are inconsistent. The fit to the new measurements by BES and BABAR 
indicates that the assumption of s = s v is not a good approximation in accordance with the goodness of fit. And 
the mixing angle is determined to be —14°, which is in good agreement with previous works. 

The content of 77 and 77' is also examined in this paper. After considering the gluonium content in the model, the 
fit to data of the world average in 2010 yields Z% = 0.002 ±0.002 and Z%, = 0.30 ±0.24, which are the contribution 
of gluonium content in 77 and 77' respectively. Although the possibility of gluonic content can not be excluded, it is 
a reasonable description for 77 in terms of pure qq meson, and no much room for a significant gluonium admixture. 
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The magnitude of gluonium contamination in rf shows that rj' has room for gluonium admixture, but the large 
uncertainty prevents us from definitely saying that gluonium content is present or not. 

As we previously stated that the latest results from BES and BABAR are not consistent with those previous 
works. The branching fractions shown in Table [T] still have large error, including statistical and systematic errors. 
The main reason is that J/ip and ip(2S) samples are not enough and the performance of detector need to be 
improved. A modern detector, BESIII[27]. has been built to meet the above requirements. Up to now, about 
2.3 x 10 8 J/tp and 1.2 x 10 8 ip(2S) events have been accumulated at BESIII, which provide a unique chance to 
study the 77-7/ mixing and further improve these measurements with much higher sensitivities. 
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